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Observation and numerical modeling of chromospheric 
evaporation during the impulsive phase of a solar flare 

S. Imada/ I. Murakami, T. Watanabe,^’^ 

ABSTRACT 

We have studied the chromospheric evaporation flow during the impulsive 
phase of the flare by using the Hinode/EIS observation and ID hydrodynamic 
numerical simulation coupled to the time-dependent ionization. The observation 
clearly shows that the strong redshift can be observed at the base of the flaring 
loop only during the impulsive phase. We performed two different numerical 
simulations to reproduce the strong downflows in FeXII and FeXV during the 
impulsive phase. By changing the thermal conduction coefficient, we carried out 
the numerical calculation of chromospheric evaporation in the thermal conduction 
dominant regime (conductivity coefficient kq = classical value) and the enthalpy 
flux dominant regime {kq = 0. lx classical value). The chromospheric evaporation 
calculation in the enthalpy flux dominant regime could reproduce the strong 
redshift at the base of the flare during the impulsive phase. This result might 
indicate that the thermal conduction can be strongly suppressed in some cases 
of flare. We also And that time-dependent ionization effect is importance to 
reproduce the strong downflows in Fe XII and Fe XV. 

Subject headings: magnetic reconnection, thermal conduction, solar flare 


1. INTRODUCTION 

Magnetic reconnection has been discussed as one of the important mechanisms for 
the rapid energy conversion of stored free magnetic energy to plasma energy. One of the 
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most famous phenomena associated with magnetic reconnection is the solar flare. A so¬ 
lar flare is a sudden brightening observed in almost all wavelengths. The total amount 
of energy released by a flare is huge and it often reaches 10^^ erg within an hour. Con¬ 
siderable efforts have been devoted toward understanding the physical mechanism of so¬ 
lar flares, and several mechanisms have been proposed over the several decades. Nowa¬ 
days, t he standard flare model which is based on magnetic r econnection, namely CSH KP 
model (jCarmichael ll964J: ISturrock Il966l: iHiravama 11974 iKopp fc Pnenman 1197611 . is 
widely believed, and modern telescope have confirmed t he predicted characte ristics from 


the model (e .g., cusp-like structure in soft X-ray images flTsuneta et ah Nl992li. reconnec¬ 


1999: 


dohva: 


Innes et ah 


2003 


1992 


ma fc Shibata 


t ion i nflows flYokovama et ah 1120011). recon nection o utflows (off limb (McKenzie &: Hudson 


Imada et ah 


1998 


Imada et ah 1120071. l2011bl) I 


Liu et ah 


20131) . on disc flHara et ah ll201l[) L plasmoid ejection 


20131) . and Coronal Mass Ejections flSvestka &: Cliver 


Recently, solar atmosphere has been focused as a space laboratory for magnetic recon¬ 
nection because of its variety in plasma condition. Actually, with the solar atmosphere, we 
can cover from low f3 (< 1) to high (3 (> 1), weakly ionized to fully ionized, and collisional 
to collisionless plasma. Observing magnetic reconnection in various plasma conditions is 
physically interesting in the held of not only solar physics but also other plasma physics. 
Especially, in the solar corona, weak Coulomb collisions can affect the plasma dynamics (e.g., 
thermal conduction, radiative cooling, ionization and recombination). In such a plasma con¬ 
dition, thermal conduction have an important role for energy transport during the explosive 
events. The plasma dynamics of magnetic reconnection with heat conduction is different 
from that wit hout heat condu ction. For example, the adiabatic slow-mode shocks predicted 


by Petschek ( iPetschek 
(iForbes et ah 


_ 1964) might be isotherm al shocks by the heat conduction effects 

19891: Yokovama fc Shibata 1997 ). Furthermore, the thermal conduction ef¬ 


fect can carry the energy to larger scale. According to the CSHKP model, the released 
magnetic energy stored in the solar corona is mainly transported t o the chromosphere by 
the thermal conduction and the non-thermal electrons flFisher 1119861) . With this process the 
chromospheric plasma is suddenly heated, and a localized high pressure region is formed. 
The pressure-gradient force drives the dense plasma to flow up toward the corona along the 
magnetic helds. This upflow is called chromospheric evaporation. 

The observational studies on evaporation have been done for several decades. A key piece 
of evidence for chromospheric evaporation is the detection of fast upflows in a high temper¬ 
ature plasma. The first observation of chromospheric evap oration was done by the Bragg 
crystal spectrometer (BCS) on P78 -1 (Feldman et ah 1980), and much work has been done 


based on the spacecraft observation (lAntonucci et ah 


I 992 I: iTeriaca et ah 120031: iMilligan &: Dennis 


2009 


I 982 I: lwatana.be 
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2013 
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as Solar Maximum Mission, Hinotori, and Yohkoh, Solar Heliospheric Observatory, Hin- 
ode, and Interface Region Imaging Spectrograph. Recently, the temperatnre dependence of 
chromospheric evaporation is intensively stndied. Inspection of the temperatnre dependent 
observation reveals several interesting characteristics of chromospheric evaporation. For ex¬ 
ample, the modern observation showed that the temperatnre of the point separating npflows 
and downflows (the flow reversal temperature) is located at ~ 1 MK. However, some obser¬ 
vations show that the flow reversal temperature locates even higher temperature during the 
impulsive phase flLi fc Ding Il201lh . Recently, to understand t hese observationa l facts the 
numerical modelings are also intensively performed. Liu et al. flLiu et ah 1 120091) discussed 
self consistently combined Fokker-Planck modeling of energetic electrons and hydrodynamic 
simulation of flare. They compared the results between conductive heating dominant case 
and direct energetic electron heating dominant case and found that the flow reversal tem¬ 
perature tend to be lower when conductive heating dominates over direct electron heating. 
Because the energy input location is different between thermal conduction and energetic 
electron precipitation, the location of the flow reversal temperature also different in height. 
They also mentioned that the unusual high temperature (2MK) downflow (~ 15 km sec“^) 
without ener getic electron is due to the thermal expansion early in the corona. Brannon 
& Longcope (iBrannon fc Longcope 1120141 ) successfully reproduce the observed flow reversal 
temperature by ID hydrodynamic simulation omitting gravity effect. Because they omit 
gravity, the solar atmosphere is not stratihed. It is not clear whether the observed flow 
reversal temperature can be reproduced when they take into account the gravity effect. The 
flow reversal temperature becomes one of the hot topics in the held of the solar physics. 


OBSERVATION 


On 6 November 2010, Hinode f Kosugi et al. 1120071) observed a solar hare (GOES M5.4, 
peak time 15:27) at the southeast (20° S, 58° E). The EUV Imaging Spectrometer (EIS) on 
board Hinode is a high sp ectral/spatial resolutio n spectrometer aimed at studying dynamic 
phenomena in the corona fjOnlhane et al. II2007I) . The hare we observed were obtained with 
hare study (HH_Flare_180xl60_v2). EIS performed a coarse raster scan (2 arcseconds slit 
and 5 arcseconds step with 8 sec exposure: sparse raster) with about 6 minutes cadence 
over this haring active region (NOAA 11121) starting before the hare, and the hare occurred 
during the scan. EIS data from the raster are processed using the EIS team provided soft¬ 
ware, which corrects for the hat held, dark current, cosmic rays, hot pixels, and slit tilt. 
For thermal reasons, there is an orbital variation of the line position causing an artihcial 
Doppler shift of ±20 km s“^ which follows a sinusoidal behavior. This orbital variation of 


the line position was corrected using the house keeping data flKamio et al. 1120101) . EIS sue- 
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cessfully observed the foot point of the flare, where the chromospheric evaporation occurs, 
during the impulsive phase of the flare. EIS obtained EUV images and line-of-sight (LOS) 
velocities estimated by Doppler shift in several emission lines (e.g., Fe XII, Fe XV, Fe XXIII). 
By using these emission lines, we can discuss the temperature dependence of the chromo¬ 
spheric evaporation upflow/condensation downflow during the flare from a few 10^ to 10^ 
K. Simultaneously, the Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly 
(AIA) instrument acquired full-Sun images with a spatial resolution of 1000 km. We us e 
the AIA 193 A passband to study the temporal evolution of the flare I Boerner et ah 1120121 ). 
Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) also observed the early 
phase of flare. The energetic electron (> 25 keV) enhancement was not observed, although 
the moderate-energy electron (~ 10 keV) enhancement was observed during the early phase 
of the flare (not shown here). The non-thermal electron acceleration might be weak in this 
event. Unfortunately, RHESSI could not observe the main phase of flare because of night 
time. 


Figure 1 shows the flare image of AIA 193 A channel at the impulsive (la) and the 
peak (lb) phase. The light curve of GOES 1.0-8.0 and 0.5-4.0 A are also shown in Figure 
Ic. The dashed vertical lines show the acquired time of the AIA images (a: impulsive, b: 
peak). During the impulsive phase, we can clearly see the two bright points at {X,Y) ~ 
(—775",—370") and (—775",—355"). It seems that these bright points represents the two 
footpoints of the flaring loop. Figure la represents the very beginning of the flare, because 
the bright points are limited in space. We can clearly observe the flare arcades at the peak 
time of the GOES light curve in Figure lb, ~4 minutes later from Figure la. There are no 
clear bright points in Figure lb any more. 

Figure 2 shows the intensity (top: a-c) and LOS velocity (bottom: d-f) images of Fe XII 
(10®'^K), Fe XV (10^'^K), and Fe XXIII (10^'^K) during the impulsive phase of the flare. The 
slit scanning starts from 15:28:08 UT, and the slit is located at the flaring region during the 
impulsive phase. We can also see the bright points in Figure 2a-b, which we already discussed 
in Figure la. On the other hand, Fe XXIII intensity image, which generally represents the 
hot flaring plasma, shows the flare loop which connects two brights in Figure 2a-b. This 
result indicates that the bright points certainly locate at the footpoints of the flaring loop. 
The LOS velocity in Fe XII at the footpoint of the flaring loop during the impulsive phase 
shows a few 10 km sec“^ downflow. The Fe XV Doppler image (Figure 2d) also shows that 
the redshift reaches 50 km sec“^. On the other hand, the Fe XXIII Doppler image (Figure 
2f) shows the ~ 50 km sec“^ upflow. We can see the temperature dependent flows of the 
chromospheric evaporation. 


From Figure 2, we hnd the dependence of flow on temperature. To understand this 
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dependence in detail, we analyze the line profiles more carefnlly. Fignre 3a-g show the 
line prohles of O VI (184.12 A, lO^-^K), Fe X (184.54 A, lO^'^K), Fe XII (192.39 A, lO^'^K), 
Fe XIV (264.79 A, lO^'^K), Fe XV (284.16 A, lO^'^K), Fe XVI(262.98 A, lO^-^K), and Fe XXIII 
(263.77 A, 10^'^K) with a single Ganssian htting, respectively. The line prohles in Fignre 
3 are obtained at the x-mark position in Fignre 2. The dashed lines show the line center 
wavelength of stationary component for each emission lines. The dotted lines represent 
the line center wavelength estimated by single ganssian htting. We can clearly see that 
the all line prohles are redshifted (dashed line < dotted line) except Fe XXIII (Figure 3g). 
Figure 4 shows the summary of temperature dependence of the chromospheric evaporation 
uphow/condensation downhow during the impulsive phase of the hare. Interestingly, the 
strongest redshift is observed in Fe XV emission line which generally represents the hot 
loops of the active region. 

Figure 5 shows the EIS observation of the peak phase of the hare. We cannot see 
the bright points anymore, and the well-developed hare arcades are clearly observed. The 
blueshift component in Fe XXIII (arrow in Figure 5c) which observed the boundary of the 
hare arcade might be related to the evaporation uphow. The strong redshift cannot be 
observed in Fe XII or Fe XV any more. Therefore, we think that the strong redshift in hot 
emission line can be observed at the footpoint of the hare loop only during the impulsive 
phase. 


3. Numerical modeling 


We have shown the Hinode/EIS observation of chromospheric evaporation uphow and 
condensation downhow during the hare. The observation shows that we can observe the 
strong downhows of hot plasma (a few MK) during the impulsive phase of the hare. In 
this section, we carry out two numerical modelings (easel: thermal conduction dominant, 
case2: enthalpy hux dominant) in order to understand the observed strong downhow during 
the impulsive phase. The calculations discussed in this paper were carried out using the 
ID version of the numerical package CANS (Coordinated Astronomical Numerical Software; 
http://www-space.eps.s.u-tokyo.ac.jp/~yokoyama/etc/cans/). 


Recently, Imada et al. f imada et ah 2nilal i]) pointed out that ionization cannot reach 
equilibrium in some part of haring region because of its fast how and rapid heating. It is 
important to take into account time-dependent ionization process when we interpret the 
observation of hare associated phenomena. Thus, we have combined the time-dependent 
ionization code with ID version of the numerical package. 
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3.1. Simulation setup 


The simulation setup is almost the same as Imada Sz Zweibel fllmada &: Zweibel II 2012 I) . 
A hxed semi circular single magnetic loop with a constant cross section is assumed. We 
assume the half-length of the loop (L) is 26 Mm. The dynamics in only half of the loop 
using a ID single fluid hydrodynamic (HD) code are calculated, because the symmetry 
about the loop top can be assumed. We used 2001 grid points in x. Grid spacing in the 
corona {x > l.Sxtr) is set to be Axj+i = 1.02Aa:j, where Xtr is the transition region height 
{xtr = 2500 km). Below the transition region {x < l.Sxtr) we use Ax = 0.01 ho, where ho 
is the pressure scale height in the chromosphere {ho = 200 km). The reflecting boundary 
conditions at x=0 and L; dp/dx = 0, dp/dx = 0, 14 = 0 are used. 

The ID hydrodynamic equation in Eulerian form are 

dp 

m 


m 


7 


P , ^ T/ 

— + jpv; 


+ 


m 

d 


dx [\7 


7 


+ A (pi4) = 0, 


( 1 ) 

d 

+P) =-P9\b 


(2) 


= -pQwVx + H - R, 

( 3 ) 

Pb rj. 

p = — pT, 
m 


( 4 ) 

go cos [(7r/2)a;/L], 


( 5 ) 


where x is the distance along a loop from the bottom of chromosphere, p is the mass density, 
V is the velocity, p is the gas pressure, T is the fluid temperature, m is the mean mass per 
particle (= 0.5mp), gi\ is the solar gravity along the loop, go is the solar surface gravity 
(2.74 X 10^ cm s-^), kB is Boltzmann’s constant, and 7 is the ratio of specihc heats for an 
ideal ga s taken to be 5/3. We use the classical conductivity for a fully ionized hydrogen 


plasma (jSpitzer Ill962l) : 


where Ko is 9.0 x 10 ^ erg s ^ K ^ cm ^ 


«|| = ( 6 ) 

The radiative loss rate {R) is expressed by 


R = p^\(p)A{T), 


( 7 ) 


where \p{p), A(T) represent the optical thickness effect on the radiative cooling efficiency 
and the radiative energy loss function, respectively. We assume Xp{p) = pd/ptsj]h.{p/pd), 
and pd = mx 10 ^^ g cm“^. As a consequence, radiative cooling is strongly suppressed below 
the transition region, where the atmosphere is optically thick. The radiative energy loss 
function, A(T), will be discussed later. 
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The flare heating is represented by the energy input rate per unit volume (H in Equation 
3.3). The flare heating is set to be symmetric and placed at the loop top and is assumed to 
be spatially Gaussian and a step function like in time: 


H{x,t) = Hoqit) 


2x/2^ 


exp 


(x — Ly 


1 + tanh 


X 


20hr 


3hn 


qit) = j 


1 + tanh 


O.lrn 


1 — tanh 


t 




O.lrn 


( 8 ) 


(9) 


where Hq, tj, Wf represent flare energy input rate (3 erg cm ^ s flare duration (240 s), 
flare region (6000 km), respectively. Tq (20 s) is the sound traveling time. 

We calculate the time evolution of iron charge states to study the effect of transient 
ionization in the chromospheric evaporation. There are many kinds of atomic species in solar 
corona, and they mainly radiate line emission in ultra-violet wavelength range by bound- 
bound proces s. Iron is the mos t dominant element for radiation in solar coronal plasma 
(a few MK) fiLandi et ah 1 120131) . Further, the recent space telescopes mainly observe the 
emission lines from iron. Therefore, we concentrate on the calculation of the time-dependent 
ionization of the iron in this study. The detail description of time-depend ent ionization 
calcula tion which is used in this study was well discussed in Imada et al. fllmada et ah 
2011 all . 


The time-dependent ionization equations for iron is 


dnf^ 


+ V . nf'v = n. + nfASa - nf‘ (af‘ + Sf')] , 


( 10 ) 


where nf ® is the number density of the zth charge state of the iron. The radiative and dielec- 
tronic recombination coefficients and the collisional ionization coefficients are represented by 
af® and Sf^, respectively. All ions and electrons have the same flow speed and temperature 
in the same location are assumed. 

We have combined ID HD code with time-dependent ionization calculation code. Be¬ 
cause we assumed that the all ions and electrons have the same temperature and velocity, 
we can calculate the time-dependent ionization (Equation 3.10) by using the temperature 
(T) and velocity (14) obtained from ID HD calculation (Equa tion 3.1-7). We cal culate the 
radiative energy loss function by CHIANTI atomic database fiLandi et ah 1120131) but with 
ionic fractions of iron calculated by Equation 3.10. Generally, bound-bound emissions are 
heavily affected by the ionic fraction. Thus, the radiative cooling may be sensitive to the 
time-dependent ionization process. We assumed all the elements except iron are in ioniza¬ 
tion equilibrium, because the dominant source of radiative energy loss is iron in the solar 
























corona. We used the usual coronal abundance flFeldman Ill992h to estimate the line emis¬ 
sions. Through the radiative energy loss function, the time-dependent ionization process can 
give feedback to the hydrodynamics. 


3.2. Simulation results 

We carry out two kinds of numerical simulations to clarify why the strong downflows can 
be observed during the impulsive phase. By changing the thermal conductivity coefficient, 
we try to study the chromospheric evaporation in the thermal conduction dominant regime 
and in the enthalpy flux dominant regime. 


3.2.1. Casel: Thermal conduction dominant regime 

First we show the result of chromospheric evaporation calculation in the thermal con¬ 
duction dominant regime. As mentioned in Section 1, the solar corona is believed to be in 
this regime. In this calculation, we have used the classical thermal conduction coefficient. 

We assume the loop is initially in hydrostatic equilibrium. The temperature along the 
coronal loop is set to be as follows; 

T(x) = r„ + i (T,.., - To) (^tanh + l) . (H) 

where Tq, T^p are the temperature of chromospheric plasma (10"^ K) and loop top (2 MK), 
respectively. Our initial condition does not satisfy thermal equilibrium. However, this is not 
sensitive to the result of our calculation, because the plasma dynamically changes by the 
energy input from the flare. We also assume the ionization equilibrium for initial condition 
everywhere. Figure 6 shows the initial condition of the calculation. The horizontal axis (x) 
shows the height along the loop from the bottom of chromosphere. The chromosphere locates 
0 < X < 2, and the corona locates 2.5 < x, respectively. From the top, plasma density (6a), 
velocity (6b), and temperature (6c) are shown, respectively. The intensity, J, of Fe XII (6d), 
Fe XV (6e), Fe XXIII (6f) in each height are also shown, respectively. The intensities are 
calculated by CHIANTI atomic database by using the electron density, temperature and nf ® 
and are normalized by the typical value (nf^/Snf® = 1 and Ue = 10^° cm“^). 

Figure 7 shows the result of 26 second later from the beginning of the simulation. From 
the top, plasma density (7a), velocity (7b), and temperature (7c) are shown, respectively. 
The normalized intensity Ipe xii (7d), Ipe xv (7f), Ire xxiii (7h) are also shown, respec¬ 
tively. The black solid line shows the result with time-dependent ionization calculation, and 
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the red dashed lines show the result with ionization equilibrium assumption. We also show 
I X Vx of Fe XII (7e), Fe XV (7g), Fe XXIII (7i) in each height, which represents the Doppler 
shift contribution parameter in each emission lines. The positive value in J x I 4 indicates 
the blueshift component, and the negative value represents the redshift component. J x I 4 
are also normalized by the typical value (I 4 = 100km sec“\ = 1, and Ue = 10^° 

cm“^). 

In Figure 7, the thermal conduction front is just reached the chromosphere, and the 
chromospheric evaporation (Vx ~ 100 km sec“^) is forming (2 < a; < 3 Mm). The pressure- 
driven downflows from the loop top, which carry a large enthalpy flux, do not reach the 
chromosphere yet. Its front locates at a; ~ 6 Mm. We can also see the condensation 
downflow (Vx ~ -10 km sec“^) at a; ~ 2 Mm. At this moment, Fe XII locates only at the 
base of the loop, and FeXV is distributed larger region (2 < a; < 5 Mm). On the other 
hand Fe XXIII is not formed yet. The ionization cannot reach equilibrium in the flaring 
region because of its fast flow and rapid heating. Note that the red dashed lines in Figure 
7 show the result of ionization equilibrium assumption, and the black solid lines show the 
result taken into account the time-dependent ionization effect. Both of the Doppler shift 
contribution parameter, I x V^, of Fe XII and Fe XV show the positive value, which means 
that the blueshift component should be observed in these emission lines. 

Figure 8 shows the result of 80 second later from the beginning. The hgure format is the 
same as Figure 7. At this moment, the chromospheric evaporation flows are well developed, 
and the highest velocity reaches 250 km sec“^. Although Fe XII locates only at the base of 
the loop, Fe XV spread much higher part of the loop (2 < a; < 7 Mm). The ionization almost 
reach equilibrium because of the high density condition of the chromospheric evaporation, 
and we can observe Fe XXIII in most of higher part of the loop. All of the Doppler shift 
contribution parameter, / x Vj,, of Fe Nil, Fe XV, and Fe XXIII also show the positive 
(blueshift) value. Therefore, the observed strong downflow in Fe XII and Fe XV emission 
line during the impulsive phase cannot be reproduced by the chromospheric evaporation 
calculation in the thermal conduction dominant regime. 


3.2.2. Case2: Enthalpy flux dominant regime 

In the previous subsection we clearly show that the chromospheric evaporation calcu¬ 
lation in the thermal conduction dominant regime cannot reproduce the observed strong 
downflow in Fe XII and Fe XV emission line during the impulsive phase. Therefore, we 
change the thermal conduction coefficient to simulate the chromospheric evaporation in the 
enthalpy flux dominant regime. In this calculation we use 10% of the classical value of the 
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thermal conduction coefficient. 

Figure 9 shows the result of 40 second later from the beginning. The hgure format 
is the same as Figure 7. The pressure-driven downflows from the loop top, which carry a 
large enthalpy flux, just reach the chromosphere. Because of the small thermal conduction 
coefficient, in this regime the enthalpy flux carry much more energy than the thermal con¬ 
duction. The highest velocity of pressure-driven downflows almost reaches 150 km sec“^. At 
this moment, Fe Xll locates only at the base of the loop, and Fe XV is distributed slightly 
higher region (2 < a; < 3 Mm). We can observe Fe XXIII in most of higher part of the 
loop. Surprisingly, in Figure 9, both of the Doppler shift contribution parameter, I x 14, 
of Fe XII and Fe XV show the negative value, which means that the redshift component 
should be observed in these emission lines. The redshift component is stronger when the 
time-dependent ionization are taken into account. 

Figure 10 shows the result of 80 second later from the beginning. The hgure format 
is the same as Figure 7. Even in this enthalpy hux dominant regime, the chromospheric 
evaporation hows are well developed, and the highest velocity also reaches 250 km sec“^. 
Although Fe XXIII distribution is slightly diherent, the result in the enthalpy hux dominant 
regime is almost the same as the result from that in the case of thermal conduction dominant 
regime. The Doppler shift contribution parameters of Fe XII and Fe XV show both positive 
and negative value, although that of Fe XXIII show the positive value. Therefore, the strong 
downhow in Fe XII and Fe XV can be observed only the beginning of the hare at the base 
of haring loop. 


4. Summary &; Discusiion 

We have studied the chromospheric evaporation how during the impulsive phase of the 
hare by using the Hinode/EIS observation. We have also preform the numerical simulation 
of ID hydrodynamics with time-dependent ionization. The observation clearly shows that 
the strong redshift can be observed at the base of the haring loop only during the impul¬ 
sive phase. We perform the two diherent numerical simulations to reproduce the strong 
downhows in Fe XII and Fe XV during the impulsive phase. By changing the thermal con¬ 
duction coefficient, we carry out the numerical calculation of chromospheric evaporation in 
the thermal conduction dominant regime (kq = classical value) and the enthalpy hux dom¬ 
inant regime (kq = O.Ixclassical value). The chromospheric evaporation calculation in the 
enthalpy hux dominant regime can only reproduce the strong redshift at the base of the hare 
during the impulsive phase. This result might indicate that the thermal conduction can be 
suppressed in some cases of the hare. We also hnd that time-dependent ionization ehect is 
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importance to reproduce the strong downflows in Fe XII and Fe XV. 

Let us discuss why the chromospheric evaporation calculation in the thermal conduction 
dominant regime cannot reproduce the strong redshift in Fe XII and Fe XV. The solar 
atmosphere are gravitationally stratifled, and the sharp temperature and density gradient 
region, so called transition region, is formed in between corona and chromosphere. In the 
thermal conduction dominant regime, the released energy by flare is conducted along the 
held line from loop top to the chromosphere. The heat from solar corona violates the 
pressure balance of the gravitationally stratifled atmosphere. The thermal conduction can 
carry only the energy, and the dense plasma is suddenly heated to form a localized high 
pressure region. The localized high pressure region should be formed at the transition region 
which has a large jump in temperature and density. This is the reason why the flow reversal 
temperature locates at the transition region temperature in the thermal conduction regime. 
On the other hand, in the case of the enthalpy flux dominant regime, the plasma energy 
and momentum are mainly transported by pressure-driven downflow from flare region to 
transition region simultaneously. Therefore, the bottom part of the flare loop is hit and 
heated by this downflow. This process can produce the strong downflow in Fe XII and 
Fe XV at the base of flare loop during the impulsive phase. 


Finally, we discuss the possibility that the thermal conduction can be suppressed in 
some cases of the flare. The classical thermal conduction in a fully ionized hydrogen plasma 
is based on the assumption that the mean free path is enough shorter than the characteristic 
temperature scale height in the solar corona. In the case that the mean free path becomes 
comparable to or even greater than the te mperature scale height , the normal diffusion ap - 
proximation for the heat flux breaks downf Cowie fc McKee 1119771: iDalton fc Balbus Ill99311 . 
Furthermore, in the case of solar flare, the large structure of magnetic reconnection can also 
contribute to suppress the thermal conduction effect. The loss-cone angle of hot electron 
in the reconnection region is small because of the large magnetic held gradient between 
upstream and downstream of reconnection region. The free streaming electron along the 
magnetic held from reconnection region to outside might be limited in a collisonless plasma. 
Suppression mechanism of thermal conduction during a solar flare is important future work. 


Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as 
a domestic partner and NASA and STFC (UK) as international partners. It is operated 
by these agencies in co-operation with ESA and NSC (Norway). This work was partially 
supported by the Grant-in-Aid for Young Scientist B (24740130), by the Grant-in-Aid for 
Scientific Research B (23340045), by the Grant-in-Aid for Scientific Research B (26287143), 
by the JSPS Gore-to-Gore Program (22001), by the JSPS Program for Advancing Strategic 
International Networks to Accelerate the Girculation of Talented Researchers under grant 













number G2602. 


REFERENCES 

Carmichael, H., in The Physics of Solar Flares, ed. W. N. Hess (NASA Special Publication 
50; Washington, DC: NASA), 451 (1964). 

Sturrock,P. A., Nature, 211, 695 (1966). 

Hirayama, T., Sol. Phys, 34, 323 (1974). 

Kopp, R. A., & Pneuman, G. W. 1976, Sol. Phys, 50, 85 

Tsuneta, S., Hara, H., Shimizu, T., Acton, L. W., Strong, K. T., Hudson, H. S., & Ogawara, 
Y. , PASJ, 44, 63 (1992). 

Yokoyama, T., Akita, K., Morimoto, T., Inoue, K., & Newmark, J., Astrophys. J., 546, L69 

( 2001 ). 

McKenzie, D. E. & Hudson, H. S., Astrophys. J., 519, L93 (1999). 

Innes, D. E., McKenzie, D. E., & Wang, T., Sol. Phys, 217, 267 (2003). 

Imada, S., Aoki, K., Hara, H., Watanabe, T., Harra, L. K., & Shimizu, T., Astrophys. J., 
776, Ill (2013). 

Hara, H., Watanabe, T., Harra, L. K., Culhane, J. L., Young, P. R., Astrophys. J., 741, 107 

( 2011 ). 

Ohyama, M. & Shibata, K. 1998, Astrophys. J., 499, 934 (1998). 

Liu, W. et al., Astrophys. J., 767, 168 (2013). 

Svestka, Z. & Cliver, E. W., History and Basic Characteristics of Eruptive Flares, in lAU 
Colloq. 133, Eruptive Solar Flares, Vol. 399, ed. Z. Svestka, B. V. Jackson, & M. E. 
Macado (New York: Springer), 1 (1992). 

Imada, S., Hara, H., Watanabe, T., Kamio, S., Asai, A., Matsuzaki, K., Harra, L. K., & 
Mariska, J. T., PASJ, 59, S793 (2007). 

Imada, S., Hara, H., Watanabe, T., Murakami, 1., Harra, L. K., Shimizu, T., & Zweibel, E. 
G., Astrophys. J., 743, 57 (2011). 



Petschek , H.E., in NASA Symposium on Physics of Solar Flares (NASA SP 50), 425 (1964). 
Forbes, T. G., Malherbe, J. M., & Priest, F. R., Sol. Phys., 120, 285 (1989). 

Yokoyama, T. & Shibata, K., Astrophys. J., 474, L61 (1997). 

Fisher, G. H., in lAU Golloq. 89, Radiation Hydrodynamics in Stars and Gompact Objects, 
ed. D. Milhalas & K.-H. Winkler (Berlin; Springer), 53 (1986). 

Feldman, U., Doschek, G. A., Kreplin, R. W., Mariska, J. T. , Astrophys. J., 241, 1175 
(1980). 

Antonucci, F., Gabriel, A. H., Acton, L. W., Leibacher, J. W., Gulhane, J. L., Rapley, G. 
G., Doyle, J. G., Machado, M. F., & Orwig, L. F., Sol. Phys, 78, 107 (1982). 

Watanabe, T., Sol. Phys, 126, 351 (1990). 

Gulhane, J. L., Fludra, A., Bentley, R. D., Doschek, G. A., Watanabe, T., Hiei, E., Lang, 
J., Garter, M. K., Mariska, J. T., Phillips, A. T., Phillips, K. J. H., Pike, G. D., & 
Alphones, G., PASJ, 44, L101(1992). 

Teriaca, L., Falchi, A., Gauzzi, G., Falciani, R., Smaldone, L. A., & Andretta, V., Astrophys. 
J., 588, 596 (2003). 

Milligan, R. O., & Dennis, B. R., Astrophys. J., 699, 968 (2009). 

Brosius J. W., Astrophys. J., 726, 133 (2013). 

Tian, H., Li, G., Reeves, K. K., Raymond, J. G., Guo, F., Liu, W., Ghen, B., & Murphy, N. 
A., Astrophys. J., 797, L14 (2014). 

Li, Y., & Ding, M. D., Astrophys. J., 727, 98 (2011). 

Liu, W., Petrosian, V.,& Mariska, J. T., Astrophys. J., 702, 1553 (2009). 

Brannon S. & Longcope, D., Astrophys. J., 792, 50 (2014). 

Kosugi, T., Matsuzaki, K., Sakao, T., Shimizu, T., Sone, Y., Tachikawa, S., Hashimoto, 
T., Minesugi, K., Ohnishi, A., Yamada, T., Tsuneta, S., Hara, H., Ichimoto, K., 
Suematsu, Y., Shimojo, M., Watanabe, T., Shimada, S., Davis, J. M., Hill, L. D., 
Owens, J. K., Title, A. M., Gulhane, J. L., Harra, L. K., Doschek, G. A., & Golub, 
L., Sol. Phys, 243, 3 (2007). 



Culhane, J. L., Harra, L. K., James, A. M., Al-Janabi, K., Bradley, L. J., Chaudry, R. A., 
Rees, K., Tandy, J. A., Thomas, P., Whillock, M. C. R., Winter, B., Doschek, G. A., 
Korendyke, C. M., Brown, C. M., Myers, S., Mariska, J., Seely, J., Lang, J., Kent, B. 
J., Shaughnessy, B. M., Young, P. R., Simnett, G. M., Gastelli, G. M., Mahmoud, S., 
Mapson-Menard, H., Probyn, B. J., Thomas, R. J., Davila, J., Dere, K., Windt, D., 
Shea, J., Hagood, R., Moye, R., Hara, H., Watanabe, T., Matsuzaki, K., Kosugi, T., 
Hansteen, V., & Wikstol, O., Sol. Phys, 243, 19 (2007). 

Kamio, S., Hara, H., Watanabe, T., Fredvik, T., Hansteen, V. H., Sol. Phys, 266, 209 (2010). 

Boerner, P., Edwards, G., Lemen, J., Rausch, A., Schrijver, G., Shine, R., Shing, L., Stern, 
R., Tarbell, T., Title, A., Wolfson, G. J., South, R., Spiller, E. Gullikson, E., McKen¬ 
zie, D., Windt, D., Golub, L., Podgorski, W., Testa, P., & Weber, M., Sol. Phys, 275, 
41 (2012). 

Imada, S., Murakami, L, Watanabe, T., Hara, H., & Shimizu, T., , Astrophys. J., 742, 70 

( 2011 ). 

Imada, S., & Zweibel, E. G., Astrophys. J., 755, 93 (2012). 

Spitzer, L., Physics of Fully lonizaed Gasses, New York. (1962). 

Landi, E., Young, P. R., Dere, K. P., Del Zanna, G., & Mason, H. E., Astrophys. J., 763, 86 
(2013). 

Feldman, U., Phys. Scr., 46,202 (1992). 

Gowie L. L. & McKee, G. F., Astrophys. J., 211, 135 (1977). 

Dalton W. W. & Balbus, S. A., Astrophys. J., 404, 625 (1993). 


This preprint was prepared with the AAS lATRX macros v5.2. 



15 


SDO/AIA 193A 
15:31:08 
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Fig. 1.— The flare image of AIA 193 A channel at the impulsive (la) and the peak (lb) 
phase. The light curve of GOES 1.0-8.0 and 0.5-4.0 A are also shown in Figure Ic. The 
dashed vertical lines show the acquired time of the AIA images (a: impulsive, b: peak). 
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Fig. 2.— The intensity (top: a-c) and LOS velocity (bottom: d-f) images of Fe XII (10^'^K), 
Fe XV (10®''^K), and Fe XXIII (10^'^K) during the impulsive phase of the flare. 
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Fig. 3.— The line profile of O VI (184.12 A, lO^-^K), Fe X (184.54 A, 10®-°K), Fe XII (192.39 
A, lO^-^K), Fe XIV (264.79 A, lO^'^K), Fe XV (284.16 A, lO^'^K), Fe XVI(262.98 A, lO^'^K), 
and Fe XXIII (263.77 A, 10^'^K) with a single Gaussian htting. 
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Fig. 4.— The temperature dependence of the chromospheric evaporation up- 

flow/condensation downflow during the impulsive phase of the flare. The LOS velocities 
are estimated by Doppler shifts. 
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Fig. 5.— The intensity (top: a-c) and LOS velocity (bottom: d-f) images of Fe XII (10^'^K), 
Fe XV (10^'^K), and Fe XXIII (10^'^K) dnring the peak phase of the flare. 
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Fig. 6.— The initial condition of the chromospheric evaporation calculation. From the top, 
plasma density (6a), velocity (6b), and temperature (6c) are shown, respectively. The inten¬ 
sity, /, of Fe XII (6d), Fe XV (6e), Fe XXIII (6f) in each height are also shown, respectively. 
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Fig. 7.— The result of 26 second later from the beginning of the simulation in the thermal 
conduction dominant regime. The normalized intensity Ipe xii (7d), Ipe xv (7f), Ire xxiii 
(7h) are also shown, respectively. The black solid line shows the result with time-dependent 
ionization calculation, and the red dashed lines show the result with ionization equilibrium 
assumption. We also show J x 14 of Fe XII (7e), Fe XV (7g), Fe XXIII (7i) in each height, 
which represents the Doppler shift contribution parameter in each emission lines. 
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Fig. 8.— The result of 80 second later from the beginning. The hgure format is the same as 
Figure 7. 
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Fig. 9.— The result of 40 second later from the beginning in the enthalpy flux dominant 
regime. The figure format is the same as Figure 7. 






































23 



X (Mm) 


Fig. 10.— The result of 80 second later from the beginning. The hgure format is the same 
as Figure 7. 







































